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Abstract. - 

We study the quantum transition from an antiferromagnet to a superconductor in a model 
for electron- and hole-doped cuprates by means of a variational cluster perturbation theory ap- 
proach. In both cases, our results suggest a tendency towards phase separation between a mixed 
antiferromagnetic-superconducting phase at low doping and a pure superconducting phase at 
larger doping. However, in the electron-doped case the energy scale for phase separation is an 
order of magnitude smaller than for hole doping. We argue that this can explain the different 
pseudogap and superconducting transition scales in hole- and electron-doped materials. 



Introduction. - High-temperature superconducting materials (HTSC) are characterized 
by strong electronic correlations which are responsible for a number of anomalous properties 
and competing phases. The occurrence of these anomalous phases is probably related to the 
close proximity of the Mott-insulating and antifcrromagnctic (AF) phases at half-filling. This 
is particularly true for hole- (p-)dopcd materials which, besides the AF and the superconduct- 
ing (SC) states, display a number of unconventional phases such as stripes [1,2], checkerboard 
structures [3], non- Fermi liquid phases, etc. It has been argued that these inhomogeneous 
phases originate from an instability of the AF phase towards phase separation [2] , which is in- 
trinsic in models with short-range interactions such as Hubbard or t — J models. On the other 
hand, long-range Coulomb repulsion enforces charge homogeneity at long distances so that 
stripes, other short-range charge inhomogeneities, or strong charge fluctuations can originate 
as a compromise between the two competing effects [2,4,5]. A number of theories support 
the notion that fluctuations of these inhomogeneities, or of related order parameters, are re- 
sponsible for the pseudogap observed in p-doped materials [2,6,7]. In this case, one would 
expect the corresponding (low) pseudogap temperature scale T* to be indirectly related to 
the phase-separation energy scale. 

The situation is not clear in electron- (n-)doped cuprates (especially in the Nd- and Pr- 
CeCuO compounds). These materials display a more continuous transition from the AF to the 
SC state [8]. Stripes have not been observed so far, although there are indications of electronic 
phase inhomogeneities [9,10]. In some mean-field calculations based upon the single-band 
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Fig. 1 - Difference between the T — grand-canonical potentials AS1 = Slsc — Qaf+sc of the SC 
and of the AF+SC phases (solid line) plotted as a function of chemical potential fj,/t. m gives the 
staggered magnetization (dashed), x the doping (dotted, scaled by a factor 4 for convenience) in the 
AF+SC phase. Data are for the Hubbard model with U — 8t and i'/i = —0.3 for hole and electron 
doping. Right: blowup of the region around the transition for the n-doped case: Here, the critical 
chemical potential u c at which Afl = is fj, c ~ 6.189, and u* ~ 6.196 (see text for its meaning). In 
both p-doped and n-doped system, m is finite at the transition point, and has a jump to m = 0. 

Hubbard model for n-doped materials, the AF gap closes continuously leading to a quantum- 
critical point (QCP) and no instabilities towards inhomogeneous phases are detected [11,12], 
suggesting that the pseudogap is directly related to fluctuations of the AF gap. On the 
other hand, calculations based upon the three-band Hubbard model yield stripes - although 
with different properties - even in n-doped cuprates [13]. 

It is an important issue to understand whether the closing of the AF gap in n-doped 
cuprates, and thus the evolution from the AF to the SC state may be related to an instability 
towards inhomogeneous phases, in qualitative analogy to p-doped materials, or whether this 
evolution is more continuous, as the phase diagram of n-doped cuprates seems to suggest [14]. 

In this paper, we address the issue of the AF to SC transition by an analysis of the 
single-band Hubbard model for parameters appropriate to n- and p-doped compounds via 
variational cluster perturbation theory (VCPT) [15-17]. VCPT is appropriate to deal with 
strongly-correlated systems, since the short-range interaction part is solved exactly within a 
small cluster. A similar analysis has been recently carried out by Senechal et al. [12]. In 
that paper, the authors show that the single-band Hubbard model is sufficient to explain the 
different shape of the phase diagram between p- and n-doped cuprates. Their results also 
suggest that the AF to SC transition is continuous and associated with a quantum critical 
point. In contrast, our results show that the issue of phase separation for the AF to SC 
transition in n-doped materials requires a high accuracy in the determination of energies, 
which is possible within the present work only thanks to two technical improvements of the 
VCPT calculations, which we are going to describe in detail below. The interplay between 
AF and SC was also studied previously via a number of numerical techniques, see, e. g. 
Ref. [18]. 

Our main results concerning the nature of the AF to SC transition (see Fig. are the 
following: In p-doped systems the situation is quite clear: the transition is first order in /x, 
i.e., there is an instability towards phase separation with a jump in doping x and staggered 
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magnetisation m at the transition (/x c in Fig. ^) . We stress again that by allowing for a more 
general spatial dependences of the order parameter, this macroscopic phase separation could 
be possibly replaced by other microscopically inhomogencous phases, such as, e.g., stripes. 
Certainly, this is expected to be the case if long-range Coulomb interaction is taken into 
account [2], On the other hand, in contrast to previous theoretical calculations, our results 
also suggest phase separation in the n-doped case, although the corresponding energy scale 
(see blowup in Fig. ^) is one order of magnitude smaller than in p-doped compounds. 

In agreement with Ref. [12], we find that the AF phase actually mixes with a weak d- 
wave SC component at small doping. A similar coexistence phase is observed, for example, in 
PrCcCuO [10], and was obtained in several calculations [19,20]. On the other hand, such a 
cohexistence obtained at the mean-field level could be just a signal for the presence of strong 
pairing fluctuations within the AF phase, as we discuss below. Upon increasing doping, there 
is a transition to a pure d-wave SC phase. At larger doping, some experiments and theories 
suggest a transition to a SC phase with different pairing symmetries in n-doped cuprates, such 
as s or mixed s and d phases [21-23]. However, this high-doping regime in beyond the scope 
of the present paper. 

Method. - We consider the two-dimensional single-band Hubbard model with nearest- 
neighbor (t) and next-nearest-ncighbor (i') hoppings, and Hubbard repulsion U. We take 
typical parameters valid for both hole- and electron-doped high-T c cuprates [24], namely 
t'/t = —0.3 and U/t = 8. Although some slight difference have been suggested between the 
parameter values (t'/t and U/t) of n- and p-doped compounds, we keep them equal in order 
to just focus on the differences due to hole- and electron- doping only. We stress that the aim 
of our calculation is a general description of the materials, so that we don't aim at carrying 
out an accurate fitting of the parameters. We have in fact checked that different values for 
the parameters (e. g., a different U/t or t'/t within 30%) or a third-neighbor hopping does 
not affect our qualitative conclusions. 

Within the CPT approach [25,26], the lattice is partitioned into disconnected clusters. 
The Hamiltonian Hq PT of the disconnected lattice can be solved by exact diagonalization. 
The interclustcr hopping Hamiltonian Vcpt is then treated pcrturbatively up to first order for 
the self-energy. VCPT (or SFA) [16, 17] generalizes this approach by decomposing H into a 
"reference" part H' and a "perturbation" V, the latter containing single-particle terms only. 
In general, one can take H' = Hq PT + AH, i. e., add to the cluster an arbitrary single-particle 
Hamiltonian AH, as long as H' is still exactly solvable numerically. In order to conserve the 
total physical Hamiltonian H = H' CPT + Vcpt, AH must then be subtracted from V [17], 
which requires V = Vcpt — AH. This extension allows for the description of symmetry- 
broken phases by introducing, via AH, a corresponding symmetry-breaking field. Notice 
that, due to this subtraction the field AH is only fictitious. However, due to the fact that the 
subtraction is only carried out in a perturbative way, results in fact do depend on AH, except 
for the case where CPT becomes exact. This apparent "arbitrariness" in the choice of AH 
is restricted by the requirement that the self-energy functional [16] (which corresponds to the 
CPT grand-canonical potential) 

n = Q' + Trln(-G C PT)-Trln(-G') (1) 

has to be stationary with respect to AH. In (JTJ, CI' and G' are the grand-canonical potential 
and the Green's function of the reference system H' , respectively, and Gcpt is the Green's 
function of the physical system H calculated perturbatively in V [17]. Physically, one can 
consider VCPT as an extension of CPT in which the interclustcr perturbation is not carried 
out with respect to the exact cluster ground state but with respect to an "optimized" state 
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parametrized by AH. This is similar in spirit to diagrammatic expansions in which the "bare" 
Green's functions are replaced by "dressed" mean-field Green's functions in some symmetry- 
broken state such as AF or SC 

Since we expect to describe both an AF and a d-wave SC phase, in our calculation 
AH contains both a staggered (h) and a nearest-neighbor d-wave pairing field (A). More 
specifically, AH = Aaf + A S c + A e , where Aaf = hJ2n( n R.,i ~ n R,i)e lQ R , and A S c = 
T Er r' v(R- ~ R')( C R,T C R',1 + h.c). Here, cr i(T destroys an electron on lattice site R with 
spin projection a. rtR j(J = c Ro .cr, i(T is the corresponding density operator, Q = (n, ir) is the 
antifcrromagnetic wave vector, and the d-wave factor r](R — R') is non vanishing for nearest- 
neighbor lattice sites only and is equal to +1 (—1) for R — R' in x (y) direction. In the SC 
term, the sum is restricted to R and R' belonging to the same cluster. 

Away from half-filling, in addition to the symmetry-breaking terms, it is necessary to add 
to the reference system a "fictitious" on-site energy, i. e. 



which plays the role of a "shift" in the cluster chemical potential with respect to the "physical" 
chemical potential fi. Without this term, which has been omitted in Ref. [12], the mean 
particle density n . = 1 =p x {x is the doping, =p corresponds to p- and n-doping, respectively) 
cannot be unambiguously determined, as different results would be obtained by evaluating 
it as n = —dil/dfi or as the usual trace over the Green's function. Of course, a consistent 
treatment of the particle density is important for an accurate analysis of phase transitions as 
a function of doping, and of phase separation. Unfortunately, the inclusion of e considerably 
complicates the numerical treatment, as now is not simply a minimum with respect to its 
parameters, but a saddle point (typically a maximum in the e direction). Physically, one can 
regard e as a Lagrange multiplier which enforces an appropriate constraint in the particle 
number [27]. Along this constraint, Q is then a minimum with respect to h and A. 

For an appropriate characterization of the phase transition, it is necessary to evaluate 
f2 with high accuracy. This is difficult to achieve by simply numerically carrying out the 
frequency integrals contained in the trace in (JJJ. It is rather convenient to exploit the fact 
that these integrals can be replaced (at zero temperature T) by a sum over the negative poles 
of GcpT(q) and of G' (see Ref. [28] for details). While the latter are directly given by the 
single-particle and single-hole excitation energies of the cluster, the former can be directly 
obtained as the eigenvalues of the matrix h a ,p(a,b) in Eq. (17) of Ref. [29] [27]. Since this 
method requires the evaluation of all cluster single-particle and single-hole states, a complete 
diagonalization of the cluster Hamiltonian must be carried out. In addition, the dimension 
of the matrix h ai p(a, b) is equal to the number of single-particle and single-hole states. This 
restricts the maximum cluster sizes that can be considered in our VCPT calculation. For this 
reason, we only consider an infinite lattice split into 2x2 clusters as a reference system. We 
stress again that this is necessary in order to achieve an accuracy that permits the resolution 
of the tiny phase-separation energy scale in the n-doped case (see blowup in Fig.^). Of course 
one could argue that finite-size corrections may introduce a even larger inaccuracy. However, 
while quantitative changes are expected for larger clusters, we observe that our result about 
phase separation is quite robust and only weakly dependent on other parameters such as t' jt 
for 0.2 < \t'/t\ < 0.4, or U/t. In addition, our calculation does show that at least such an 
accuracy is required in order to address the issue of phase separation in the n-doped case. 

Results. - Our calculation, thus, proceeds as follows. For a given value of the variational 
parameters e, h, and A we carry out a complete diagonalization of the cluster Hamiltonian (due 
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to the presence of the pairing field, particle number is not conserved), evaluate il as described 
above, and look for a stationary solution as a function of the variational parameters. Quite 
generally, different stationary solutions can be found, corresponding, for example, to different 
phases. In this case, the minimum fl selects the most stable phase. Near half-filling, the 
two most stable solutions are a coherently mixed AF+SC and a pure SC phase, in agreement 
with Ref. [12]. Results for these solutions are plotted in Fig. for the p- and n-doped cases. 
The AF order parameter m is plotted as a dashed line for both n- and p-doped systems as 
a function of /i. m monotonously decreases upon going away from half-filling. In Fig. ^ 
we also plot the doping x and the difference between the T = grand-canonical potentials 
Sfl = &sc — ^af+sc °f the SC and the AF+SC solutions. While at lower doping the AF+SC 
phase is stable, a crossing to the SC phase occurs at a critical chemical potential /i c for which 
m is still nonzero, i. c. the transition is first order as a function of fj, both for the p- as well as 
for the n-doped case, as anticipated. Although, from Fig.^ this effect appears to be very small 
for the n-doped case, we stress that our calculation was carried out with sufficient accuracy 
to resolve this energy difference. At this value of //, the doping x c of the AF+SC phase is 
different from the one x' c of the pure SC phase so that there is a jump Ax = x' c — x c > in 
the doping at [i c indicating phase separation between a weakly doped AF+SC and a higher 
doped SC phase. Eventually, for some /i = fj,* (further away from half filling than /i c ) the 
AF+SC solution ceases to exist. In principle, an unstable AF+SC solution continues to exist 
for x > x c . However, in this region /i is a decreasing function of the particle density, which is 
an equivalent indication of phase separation. 

The difference A[i = \fi* — fx c \ can be considered as a measure of the characteristic phase- 
separation energy per particle, i.e., it is proportional to the energy barrier between the two 
doping values x c and x' c . As one can see from Fig. ^ m the p-doped case this energy scale 
(Afi/t W 5 • 1CP 2 ) is about an order of magnitude larger than in the n-doped case (A/x/t w 
7 • 1CP 3 ). Taking typical values for the energy unity t (t ~ 0.25eV), this corresponds to a 
temperature scale of ~ 125A" in p-doped, and ~ 18AT in electron- doped cuprates, which is 
roughly of the order of the corresponding (low-energy) pseudogap temperature scales observed 
in these materials [30]. Also the discontinuity Ax (not shown) is larger in the p-doped case 
(Ace = 0.11, vs. Acc = 0.08 for n-doped), although the difference is not significant. At 
the same time, the doping at which the AF+SC solution is destroyed is much larger in n- 
doped (x c « 0.13%) than in p-doped systems (x c ~ 0.03%), in qualitative agreement with 
experiments [8,14]. Therefore, our results suggest that the AF+SC to SC transition as a 
function of /i is clearly first order in the p-doped case while weakly first order for n-doping. 

Single-particle spectrum. - The different behavior between n- and p-doped systems can 
be understood by considering the doping evolution of the single-particle spectrum from the 
AF insulator to a SC displayed in Fig. [3 In agreement with experiments [31,32], in p-doped 
cuprates doped holes first enter at (tt/2,tt/2) [33]. On the other hand, in n-doped materials 
doped electrons initially form pockets around (7T,0) [11,31,32]. At (tt, 0)the density of states 
is larger, and the fact that \i lies within the SC gap apparently stabilizes the AF solution for 
a larger doping range, allowing for the AF gap to decrease more gradually in the n-doped 
case. The transition to a non-magnetic solution appears to occur as soon as the Fermi energy 
touches the top of the band at (tt/2,tt/2). At this nodal point the SC gap is zero, so that 
doping into nodal particles apparently destabilizes the AF solution. This appears to be the 
reason why in the p-doped case the AF solution, where hole are first doped near (tt/2, 7r/2)(scc 
also Ref. [11]), is stable only in a smaller doping range. Numerically, we observe the following 
behavior: as long as \i remains below (above for p-doped) the band around (tt/2, tt/2), there is 
only an absolute minimum of fl at a finite value h op t of the staggered field h (we keep the other 
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Fig. 2 - Single-particle spectrum for the hole-doped (a), and for the electron-doped system (b) just 
before the AF+SC to SC transition, i.e., for dopings x ~ 0.03 and x ~ 0.13, respectively. 



two fields A and e at their saddle point). As soon as fj, enters the band at (tt/2, it/2), fl(h) also 
develops a local minimum at h = 0, which becomes rapidly lower in energy than the minimum 
h = h opt . Eventually, the local maximum lying between the two minima at h = and h = h op t 
merges into the minimum at h — h opt so that this latter minimum disappears. The observation 
that doping into (ir/2, 7r/2)makes the AF phase unstable suggests that the occurrence of phase 
separation is generic and quite independent on specific parameters (unless, of course, one adds 
some doping-dependent potential). Finally, let us mention that the Fermi-surface evolution 
of the n-doped system as a function of x (not shown) qualitatively reproduces the ARPES 
experiments [32]: For small doping, we obtain electron pockets around (it, 0), while for larger 
doping our results display an evolution to a large Fermi surface centered around (tt, tt). 

Discussion and conclusions. - The VCPT method exactly treats fluctuations up to the 
range of the cluster size, so that the question arises whether the SC solution we (and also 
others [12, 19,22]) obtain within the AF phase is a true long-range SC phase or whether it is 
only a signal of strong pairing fluctuations within the AF phase leading to a SC pseudogap. 
The latter hypothesis could be supported by the fact that results obtained with different 
cluster sizes [12]seem to indicate a size dependence of the SC order parameter, and by the 
fact that the SC order parameter is about a factor three smaller in the AF+SC phase than in 
the pure SC one. The presence or not of such microscopic cohexistence phase may depend on 
material details. Certainly, our results suggest that the SC gap (or pseudogap in the case of 
fluctuations) A is important in order to stabilize the AF phase in n-doped materials. On the 
other hand, A hardly has an effect in p-doped systems since it does not produce any gap at 
the Fermi Surface, as doped holes first enter near the nodal points. 

In conclusion, our calculations suggest that the destruction of the AF phase in n-doped 
systems is associated with a tendency towards phase separation, similarly to p-doped materi- 
als, although with an energy scale which is one order of magnitude smaller. In the presence 
of Coulomb interaction, this tendency is expected to produce a microscopically segregated 
phase, such as stripes or similar inhomogeneous structures [13] Following the idea that the 
pseudogap temperature is related to the formation of such segregated phases [2], our results 
provide reasonable orders of magnitudes for the pseudogap temperature scale, and, in par- 
ticular, provide a possible explanation why this is much smaller in n-doped than in p-doped 
cupratcs. The results of the present paper support the idea that the asymmetry between hole- 
and electron-doped materials can be understood within a simple single-band Hubbard model, 
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at least concerning the portion of phase diagram we studied here. Of course, it remains an 
open issue, whether a three-band model may be necessary for a complete description. In 
fact, alternative theories propose that Ce- doping in n-dopcd compounds, such as NdCeCuO, 
occurs in a total different way. The idea there is that doped Ce atoms do not add electrons 
to the Cu02-layers, but rather are responsible for a quenching of Cu moments [34] . 

* * * 
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